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We investigate the relation between mobility and diffusivity for Brownian particles under steady 
shear near the glass transition, using mode coupling approximations. For the two directions perpen- 
dicular to the shear direction, the particle motion is diffusive at long times and the mobility reaches a 
finite constant. Nevertheless, the Einstein relation holds only for the short-time in-cage motion and is 
violated for long times. In order to get the relation between diffusivity and mobility, we perform the 
limit of small wavevector for the relations derived previously [Phys. Rev. Lett. 102 (2009), 135701], 
without further approximation. We find good agreement to simulation results. Furthermore, we split the 
extra term in the mobility in an exact way into three terms. Two of them are expressed in terms of mean 
squared displacements. The third is given in terms of the (less handy) force-force correlation function. 

§1. Introduction 

In thermal equilibrium, the Einstein relation 1 ' for a Brownian particle (a colloid) is 
the most known application of the fluctuation dissipation theorem (FDT). It connects the 
equilibrium mobility ^ e \t) of the particle to its equilibrium diffusivity D^ e \t), 

^\t)t = k B T^\t) 1 (1-1) 

with thermal energy ksT. The diffusivity on the left hand side of this equation is defined 
via the mean squared displacement (MSD) of the particle, 

2Z?( e )(t)t = ((s(t) - s(0)) 2 ) = 5s 2 (t). (1-2) 

s denotes an arbitrary coordinate of the position of the tagged particle. 5s 2 (t) denotes the 
1-dimensional MSD of the 3-dimensional system in equilibrium, which is independent of 
direction because of isotropy. (...) denotes an ensemble average, which is defined in detail 
below. The mobility on the right hand side of Eq. (11-11) is probed by a constant test force 
F which stalls to act on the particle at t = 0, and defined as the ratio of the mean velocity 
and the force, see also below (Eq. ( 12-91 )). For a single colloid in a solvent, the MSD grows 
linearly in time and the mobility is time independent, on the Brownian timescale considered 
here. 2) For higher densities of the colloids, the MSD shows the characteristic two-timescales 
scenario with short and long time diffusivities 2) ' ' and both diffusivity and mobility are time 
dependent. At even higher densities, colloidal dispersions exhibit slow cooperative dynam- 
ics and form glasses. In these, long time diffusivity and mobility vanish since the particles 
are trapped by the surrounding particles. 4) These metastable soft solids can easily be driven 
into stationary states far from equilibrium by already modest flow rates. Under shear, the 
system recovers ergodicity 5 ^' 6) and the particles recover finite diffusivities 7 ^ 11 ' and, as will 
be shown, finite mobilities for the directions perpendicular to shear. The MSD in shear di- 
rection in contrast grows cubically in time. 7) ' 12) ' 13) In this paper, we will focus on the two 
directions perpendicular to the shear direction. Although these show diffusive behavior, the 
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relation CED does not hold any more in general and much work is devoted to understanding 
the relation between diffusivity and mobility in non-equilibrium. The violation of the FDT 
is generally described by the fluctuation dissipation ratio (FDR) Xf(t). It depends on the 
considered observable / and for the case of the Einstein relation, we define the FDR via 

X il {t)^D(t)t = k B Tn{t). (1-3) 

The equilibrium-FDT in Eq. (1 1 - 1 b is recovered with (t) = 1. In non-equilibrium, X^(t) 
deviates from unity. This is related to the existence of non-vanishing probability currents. 
FDRs are hence considered a possibility to quantify the currents and the deviation from 
equilibrium. 14 ) 

The FDT-violation of colloidal glasses in a non-equilibrium steady state under shear 
for auto-correlation functions was discussed in previous papers. 13 )' 15 ) The general finding 
is that at short times, the equilibrium-FDT holds with Xf = 1, whereas at long times, 
the equilibrium-FDT is violated in a special way: the FDR lim^ 00 Xf(t) = Xf(r) = 
Xf is time-independent for the whole final relation process of the driven glass, r is time 
rescaled with the timescale of the shearing. This finding is in agreement with spin-glass 
predictions 16 ^ as well as detailed simulations. 9 ) While in Ref. 9), the value of X{ was found 
to be independent of observable /, which lead to the notion of an effective temperature T e g, 
we indeed found a slight dependence of Xf on observable. 13 )' ) 

In the case of the Einstein relation, the MSD grows linearly in time at long times and 
the mobility is constant, as we will see. Due to this, the FDR naturally approaches a constant 
value (different from unity) at long times. It is still interesting to see whether there is a sharp 
transition in the parametric plot of the two quantities from short to long time behavior as 
was found for the correlators 9 )' 15 ) and in simulations also for the Einstein relation. 9 ) The 
Einstein relation for sheared glassy systems has also been studied in Ref. 17). 

The Einstein relation is probably the FDT-example which is most easily measurable in 
experiments. Its non-equilibrium version has indeed been investigated for a single driven 
Brownian particle. 18 ) As far as known to us, the FDT-violation for the sheared system has 
not been studied in experiments before. 

In this paper, we present the Einstein relation for colloidal suspensions under shear. 
We will therefore perform the limit of small wavevector for the relations presented previ- 
ously, 13 )' 15 ) without further approximation. In Sec.|2j we give the introduction to the system 
under consideration and derive exact expressions for the quantities of interest. In Sec. 
these exact expressions will be made calculable according to previous approximations. In 
Sec. [4] we show and discuss our final results and we close with a summary in Sec. [5] 

§2. Microscopic starting point 

We consider the same system as in Ref. 15) and give its description here for com- 
pleteness. It contains N spherical Brownian particles of diameter d, dispersed in a solvent 
in volume V. The particles have bare diffusion constants Dq = ksT^iQ, with mobility 
fiQ. The interparticle force acting on particle i (i = 1 ... N) at position is given by 
Fj = — d / dviU{{rj}), where U is the total potential energy. We neglect hydrodynamic 
interactions to keep the description as simple as possible. These are also absent in the com- 
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puter simulations 9 ) to which we will compare our results. 

The external driving, viz. the shear, acts on the particles via the solvent flow velocity 
v(r) = 7yx, i.e., the flow points in x-direction and varies in y-direction. 7 is the shear rate. 
The particle distribution function \P(r = {r-j}, t) obeys the Smoluchowski equation, 2 )' 5 ) 

W{r,t) = fi9{r,t), 

Q = Q e + 5Q = 9i ■ [di - Pj - k • n] , (2-1) 

i 

with k = 7xy for the case of simple shear, fl is called the Smoluchowski operator and it is 
built up by its equilibrium part, f2 e = ^ di ■ [di — Fj] of the system without shear and the 
shear term 5f2 = — ^ di ■ k • r.;. We introduced dimensionless units d = ksT = Do = 
1. There are two time-independent distributions, the equilibrium distribution \P e oc e~ u , 
i.e., Qe^e = and the stationary distribution & s with = 0. Ensemble averages in 
equilibrium and in the stationary state are denoted 

(...}= / dr* e {r) . . . , (2-2a) 



drv s (r)..., (2-2b) 



respectively. In the stationary state, the distribution function is constant but the system is 
not in thermal equilibrium due to the non-vanishing probability current, which gives rise to 
the violation of the equilibrium-FDT. 5) ' 15) 

2.1. Mean squared displacement 

While the coherent, i.e., collective dynamics of shear melted glasses has been discussed 
in detail, )' 19) we focus here on the MSD of a tagged particle. Its general properties under 
shear for the different directions will be presented elsewhere. 7 )' 8) The MSD follows from 
the q — > limit of the incoherent density correlation function. 4 )' 20 ) In the following, we con- 
sider the (1 -dimensional) MSD parallel to the unit-vector ej_, one of the two unit-vectors 
and spanning the plane perpendicular to the shear direction. The long time diffusivity 
(and also the mobility) is slightly anisotropic in this plane, 10 )' 11 ) but the relation between 
diffusivity and mobility derived below will be independent of direction. It couples the dif- 
fusivity of a certain direction to the mobility of the same direction. 

From the different time dependent correlation functions 6 )' 21 ) after switch-on of steady 
shear, we distinguish different MSDs. They will enter the final formula for the stationary 
mobility and we introduce them briefly. In the stationary state, one measures the stationary 
MSD 5rf(t) as the q — > limit of the stationary tagged particle density correlator C q (t) = 

^ e -tq.r. e l3ti e iq(t)*^(7) foj . ^ direction q = ^ 

2 - 2(7 i (t) 

2D(t)t = 5rf(t) = lim ^ — . (2-3) 

<?-+o q l 

r s is the position of the particle and q(t) = q — q • nt is the advected wavevector which en- 
ters through translational invariance of the considered infinite system. 5 *' 6 ) For the directions 
perpendicular to shear considered here, the wavevector is time-independent. Stationary dif- 
fusivity D(t) as well as stationary mobility fi(t) should carry an index i which is suppressed. 
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If the MSD measurement is started a (not too large) period t w after switch-on of steady 
shear, one measures the two-time MSD Srf(t,t w ), where t is still the correlation time, 
i.e., 5rf(0,t w ) = 0. For the special case of t w = 0, all quantities are denoted tran- 
sient. The transient MSD follows from the transient density correlator C q (t, t w = 0) = 

^ e -iq-r 3e l7 1 'i e jq(i)-r s \ as 

2-2C < (t,0) 

<5r? (t , 0) = lim . (2-4) 

g->0 q z 

For finite we use the integration through transients (ITT) approach 5 ) in order to express 
the solution of Eq. (12- lb a time t w after switch-on as 

&(t w ) = & e + / dse ns f2V e . (2-5) 



When performing averages with \P(t w ), one uses partial integrations to let the operators in 
(I2-5I ) act on what is averaged with \P{t w ). The ITT approach has proven useful for deriving 
the stress under time-dependent flow. 22 '' 23 ^ The two-time MSD 6r?(t, t w ) follows, 



r 1l \ T (i, t % 

br i [t, t w ) = lim f 

9^0 q z 



lim 

g->0 



(2-6) 



<?xy = — FfVi is a stress tensor element which followed from Q^ e = 7cr. r j / <P r e . Oper- 
ators act on everything to the right, except for when marked differently by bracketing. For 
very long waiting times, one has \P(t w — > oo) — > ^ in Eq. (12-51) . and Srf(t, oo) = 8rf(t) 
holds. In the quiescent system, i.e., without shear, one measures the equilibrium MSD, 
which follows from the equilibrium correlator Cq S \t) = (e~ tq ' rs e flet e !q ' rs }, 

5s 2 (t) = ]im 2 ~ 2C j e)it) . (2-7) 
g^o q z 

The un-sheared system is isotropic and only the magnitude of the wavevector enters in 

c [ q e \t). 

2.2. Mobility 

Let us derive the formally exact expression for the stationary mobility /j,(t) of the tagged 
particle in direction e l ± . Therefore we consider the susceptibility for tagged particle density 

fluctuations in this direction, x qe l ± (*) = ( 8e Q r ^ - ■die Q ^ t e iqe ''^' ra )^\ as found by Agarwal 
in 1972, 24) see Ref. 5) for a derivation. The mobility fj,(t) for direction ej_, we are seeking 
here, follows from the q — ► 0-limit of the susceptibility above, 

Xae 1 (*) 

H(t) = lim -^i— . (2-8) 
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Physically, the mobility fi(t) is the ratio of the average velocity v p = v p e l ± of the tagged 
particle under the external force F(t) = Fe l ± 0(t) with step function 0, 



F 



(2-9) 



(. . . denotes an average in the steady state which has been perturbed by the force F. 

For the sheared system, the mobility is always finite in contrast to un-sheared glasses. 

In Refs. 13), 15), we presented the exact splitting of the susceptibility into four terms. 
The first represents the equilibrium-FDT and the extra term is split up into three terms. The 
limit of q — > can be done in a straight forward manner to yield the mobility (compare 
Eq. (20) inRef. 13)), 



=o 



7 

lim — ~ 



ds(a xy e 



fit; 



(fill 



-ige*,-r s \ p QU Jqe\ r 



(2-10) 



Again, the first term on the right hand side of Eq. d2T0| ) is the equilibrium Einstein relation 
for 1 -dimensional diffusion, the other three terms hence correspond to the violation of the 
equilibrium Einstein relation. We have identified all but one terms of the susceptibility under 
shear with measurable mean squared displacements. The contributions of the different terms 
in x q (i) were additionally estimated with full microscopic mode coupling approximations 
inRef. 13). 

The last term, which has yet no clear physical meaning for finite q, can be connected to 
the force-force correlation function for q — ► 0: Performing the limit of small q, many terms 
vanish due to J °° ds(a xy e n s l) = and we are left with 
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g^o 2q 2 J 



ds(o- xy e aU (&e 
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ds{a xy e 



Fle nU ri 



iqe' ± -r s ^ e QU e iqe l ± -r s ^ 



F l r l 



with Ft = F s ■ e l , , rl = r s ■ e\ and 



d_ 



ds{a xy e QU Fie nH Fl). 



(2-11) 



(2-12) 



Eq. (12421 follows with Qh^ = F*. Thus, /i 4 (t) can be expressed in terms of the force-force 
correlation function 



C Fl (t,t w ) = {Fie QH Fl)+^ I ds(a xy e^ s Fle^Fl) 

Jo 



(2-13) 



in the following way 



H 4 (t) = F(t) - J^(0) = / dt' 

Jo 



7 



ds{a xy e Q ' s Fie QH 'Fl) 
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= - dt' [C Fi {t',0) - C F i(t',oo)] . (2-14) 
^ Jo 

The last term in Eq. (12-101 ) is hence connected to the difference of transient and stationary 
force-force correlation function. For long times, the force correlation decays and (14(f) 
reaches a constant. All terms of the mobility are now connected to well defined correlation 
functions in an exact way. The first three terms in Eq. (12- 10b yet have a different quality 
compared to (ii(t), which we want to stress concerning the discussion about different forms 
of non-equilibrium FDTs: 25 )~ 27) The mean squared displacements are much easier to be 
determined experimentally or in simulations than the force correlation function, which can 
only be found when the particle positions can be resolved accurately. We note that the 
stationary force correlation function in Eq. (12-141) is equal to the one in Eq. (13) of Ref. 28). 
Performing the limit of q — > starting from Eq. (18) of Ref. 13), we get directly 

m = iirt 6r * {t) ~\t dt ' c ^ ti ' °°) + \ (2 - i5) 

With this, Eq. (13) of Ref. 28) is reproduced. We want to emphasize again that we judge 
Eq. (12-101 ) more useful compared to Eq. (12-151) . 

§3. Approximations 

In this section, we will derive approximations in order to find closed expressions for the 
exact relations above. First, we treat the mean squared displacements, and then the mobility. 

3.1. Mean squared Displacements 

3.1.1. Equation of motion for the transient MSD 

Within MCT-ITT, the mode coupling approach for sheared suspensions, the general 
strategy is to derive equations of motion for the transient quantities. These have the advan- 
tage that they are the input for the generalized Green-Kubo relations, 5 )' 29 ) e.g. for the shear 
stress. 6) The transient quantities are also more handy to be analyzed since they show e.g. the 
same plateau values as the corresponding functions in equilibrium. 6 )' 19 ) In a second step, 
two-time and stationary quantities are derived via the ITT formula, Eq. (12-51) . The separation 
parameter e describes the distance from the glass transition density. It is positive for glassy 
and negative for fluid states. The derivation of the equation of motion for the coherent (i.e. 
collective) transient correlator C™ h (t,0) = (p\e QH p q(t) ) / '(p*p q ) , with p q = J2i eiq ' r * 
has been presented in Ref. 6). Ref. 7) presents the incoherent, i.e., single particle dynam- 
ics, which will also be published in a forthcoming paper. 8 ) The equation of motion for the 
transient mean squared displacement for direction e\_ is similar to the one for the quiescent 
MSD 4) ' 30) and reads 7 ) 

Srf(t,0)+ f m (j,t-t')5rf(t,0)=2t. (3-1) 
J o 

While Eq. (13-11 ) is still exact, one has to make approximations in order to find the memory 
function 7770(7, i). The MCT-ITT route leading to its numerical evaluation is presented in 
Appendix A. Then, 5rf(t, 0) and 5s 2 (t) (with mo(0, £)) are derived with Eq. (13-11 ) . 7710(7, t) 
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and a below are the two quantities needed from the Appendix, all other quantities follow 
from the equations presented in the main text. See Ref. 31) for a review on MCT-ITT. 

3.1.2. Two-time MSD 

Having found the transient MSD, we now derive approximations in order to go from 
the transient to the two-time MSD. In Refs. 7), 15), 21), this approximation was presented 
for the two-time correlator. The corresponding approximation for the MSD follows by the 
limit of q —> 0. The result, whose derivation will be presented in Ref. 21), reads 



(3-2) 



7 9t v " 

The factorization of t w - and t-dependent terms yielded the pre-factor 

rt w . 

H t w)=j ds(a xy e Q s a xy )/{a xy a xy ). (3-3) 
Jo 

It contains in the numerator the shear stress a(t w ) = 7 J^ w ds(a xy e n ^ s o~ xy ) measured after 
switch-on. It grows with t w 32 ' ) and, for t w — > 00, reaches the familiar steady shear stress 
measured in 'flow curves' as function of shear rate. According to Eq. (13 -2b . the two-time 
MSD equals the stationary one once the stress has reached its steady value. 

Eq. (13-21 ) is exact for small waiting times, where a(t w ) = jt w + O(i^) holds and 
Eq. (13-21 ) contains the first two terms of the Taylor expansion in t w . For longer waiting 
times, it holds also quite well. 21) The important connection of the waiting time derivative in 
Eq. (13-21 ) to easier accessible time derivatives is given by ^' -* 



l(Srf(t,0)-Ss 2 (t)). (3-4) 



Eq. (13-41 ) has been successfully tested quantitatively. 21 ) Putting Eqs. (13-21 ) and (13-41 ) together, 
we finally have 

5rf(t,t w ) « 5rHt,0) + (<5r 2 (t,0) - 5s 2 (t)) . (3-5) 

In this relation, we still need to know the normalized shear stress a(t w ), which in MCT-ITT 
is expressed in terms of the coherent density correlator. Only a = a{t w — > 00) is needed for 
the graphs below, see Appendix B for its derivation. From Eqs. (13-11 ) and (13-51 ). the diffusive 
long time limit of the MSD follows, 

lim 5r 2 (t, t w ) = 2Dt = — ^ (£>0 ^ 0) ( 3 . 6 ) 

7710(7, 2 = 0) is the Laplace-transform £ {mo(7, t)} (z) = J °° dte~ zt mo(j, t) at z = 0. It 
scales in the liquid with the a-relaxation time r a and with | -y | 1 in glassy states. In the glass 
(e > 0) the scaling with shear rate in Eq. (13-61 ) follows. The long time limit of 5rf(t, t w ) 
is ^-independent, since the second term in Eq. (13-51 ) approaches a constant. This must of 
course hold because at long times, the two-time MSD has to follow the steady diffusivity D 
for all L„. 
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Fig. 1. Stationary MSD and twice the integrated mobility fi'(i) = / dt'(i(t') in direction eY for a glassy state 



(e = 10~ d ). Shear rates are j/T = 10"" 



from right to left. 



3.2. Mobility 

We want to derive approximations for the exact Eq. (12-101) following Refs. 13), 15) for 
general observables /. The final non-equilibrium fluctuation dissipation relation derived 
there compares well to the simulation results from Ref. 9). In mode coupling estimates, it 
was found that the last two terms in Eq. (12-101 ) have different sign and almost cancel each 
other. 13) We hence ignore them here as was done in Refs. 13), 15). A key approximation is 
again the connection of the waiting time derivative to time derivatives, Eq. (13-41) . We finally 
have the approximate expression for the mobility in terms of mean squared displacements, 

* lP rHt) ~ ^ (< ^ M) " ^ • (3 ' 7b) 
This is the main result of this paper. 



1 

4 dt v 

id_ 

Idi 



(3-7a) 



54. Results 



4. 1 . Mobility at short and long times 

We are now ready to discuss our results. First, we want to illustrate the outcome of 
our Einstein relation Eq. (|3-7b| ) using the equations of Sec. 13.11 as input. The transient 
MSD is given as solution of Eq. (13-11 ). the stationary MSD follows then with Eq. (13-51 ) and 
t w —s- oo. The mobility is calculated with Eq. d3-7b| ). We present the time integrated version 
of Eq. (!3-7bb for convenience. Fig. Q] shows the stationary MSD together with the time 
integrated mobility for a glassy state at different shear rates. For short times, i yb 1, we 
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Fig. 2. Parametric plot of integrated mobility fx'{t) = / * dt'fi(t') versus stationary MSD for a glassy state 
(e = 10 -3 ). The dashed dotted line shows the 'ideal' situation, where Srf(t, 0) ~ Sr^it) is used in 
Eq. d3-7bt . Inset shows the FDR as function of strain jt for the same susceptibilities. 



see that the equilibrium Einstein relation holds. We have 5r 2 (t, 0) = 5s 2 (t) + 0(7i) 7) ' 19) 
and with Eq. dTTbb . 

tit) = \§- t Sr i(t) Mt) = 1 ¥ < i- (4-D 

F° r t > |t| » the mobility is smaller than expected from the Einstein relation. The para- 
metric plot in Fig. |2] shows a rather sharp transition from short time to long time behavior 
with a straight line at long times corresponding to a constant FDR X^. In the inset of Fig. 12 
we see that the transition from the short to the long time value of nevertheless takes two 
decades in time and the strain, at which this transition happens depends strongly on shear 
rate. This is not apparent in the parametric plot. In Fig. [2] we also show the mobility from 
Eq. (13 -7bl > with transient replaced by stationary MSD, from which a sharp kink in the para- 
metric plot and a sharper transition in the inset follows. This approximation was referred to 
as 'ideal X = i-law' in Ref. 13). All our findings are reminiscent of the FDT-discussion for 
finite q. For t > | | 1 with 7 — > 0, we have Ss 2 (t) = const, in Eq. (13-7bl) and the mobility 
is given by 

lim M (t) = ~ (5r 2 (t) - l -5r 2 (t^ t > (4-2) 

For very long times, t S> It 1 )" 1 , the mobility in Eq. (14-21 ) reaches the constant £l oc I7I 
(glass). The proportionality to shear rate follows with Eq. (13-6I ). The long time Einstein 
relation under shear connecting D and (1 is then given by (we restore physical units only for 
this equation), 

X^t) = k B Tfi, (4-3) 
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2 4 6 8 

log 10 rt 



Fig. 3. Stationary MSD and integrated mobility fj! (t) — f* dt'fi(t') for a fluid state (e = — 1CP 3 ). Shearrates 
are j/T = io~ 8, ~ 6 '~ 4 '~ 2 from right to left. Inset shows the parametric plot for the different shear rates. 
For the smallest shear rate, it is almost indistinguishable from the equilibrium-FDT line. 



with 



2 dt 



(«5rf(i)-|*rf(t,0)) _ i 



Xu = lim ^ Va " ' = n" ( 44 ) 



t- 



2 dt 



The last equality followed from the equality of transient and stationary MSDs for jt S> 1. 

A comment concerning the nontrivial appearance of X^ = ^ for all glassy states in 
Eq. (14-41) is in place. Let us briefly recall the approximations which lead to Eq. (14-41) . First, 
the second term in the mobility, Eq. ( 12-101 ), the waiting time derivative, is expressed in 
terms of time derivatives, see Eq. (13 -Vbb - This relation has been tested in switch-on simu- 
lations 2 ^ and holds quantitatively for two different super-cooled liquids, at least for small 
shear rates. Second, the third term in Eq. ( 12-101 ) is neglected. It is the difference between 
the time derivatives of transient and stationary MSDs. As stated above, at long times, the 
MSDs must follow D independent of waiting time, and the third term vanishes making this 
approximation exact. Third, the last term in Eq. (12- 10b . ^^{t), is neglected. To our current 
knowledge, it is the only term which can give a deviation from the value of ^ at long times, 
its limiting value is given with t — > oo in Eq. (|2- 14b . This has to be addressed in the future. 
For finite q, we found a correction to the small shear rate glass-value X = i. 13 )> 15 ) The 
correction followed from the difference between transient and stationary correlators, this 
difference is absent for the MSDs at long times. As noted before, many spin models yield 
X = \ at the critical temperature. 33 ' 

Next, we consider the FDT-violation for the MSDs on the fluid side. We find the ex- 
pected behavior, i.e., for small shear rates, the equilibrium-FDT is restored for all times, see 
Fig. [3] All our findings are in agreement with the simulation results in Ref. 9). 
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Fig. 4. Comparison to simulation data for stationary MSD and integrated mobility = f dt'fj,(t') in 

the neutral direction at temperature T — 0.3 (T c — 0.435) and 7 = 10~ 3 . Circles and squares are the 
simulation data (including units) from Fig. 16 in Ref. 9). The full line is the integrated mobility calcu- 
lated via Eq. B-7bb . The dashed line shows the integrated mobility from Eq. Q-7bt with approximation 
5z 2 (t, 0) « Sz 2 (t). Inset shows the different MSDs, see main text. 



4.2. Comparison to simulations 

Following the discussion of the approximations for the different terms, it is interesting 
to compare Eq. (!3-7bb directly to the simulation results in Ref. 9). This is done in Fig. [4] We 
need transient as well as quiescent MSDs from the simulations as input which are not avail- 
able. Therefore, we construct a quiescent MSD which is constant for long times starting on 
the plateau of the stationary MSD. Approximating transient and stationary MSD to be equal, 
we then get the dashed curve for the mobility in Fig. 01 without adjustable parameter. In a 
second step, we can take the difference between stationary and transient MSDs into account 
with Eq. ( 13-51 ). We use the same value a = 0.01 as in Fig. 4 in Ref. 15). There, we compared 
the susceptibility for density fluctuations at wavevector q = 7.47 to the simulation data of 
Ref. 9). Since the simulation data in Fig. [4] and Fig. 4 in Ref. 15) are for the same system 
at equal shear rate and temperature, the value of the fit-parameter a should be the same in 
both figures. The constructed quiescent and the calculated transient MSDs are shown in the 
inset. The resulting mobility fits very well to the simulation results. It is interesting to note 
that the long time mobility (jt S> 1) as calculated from d3-7bb is independent of our choice 
of the height of the non-ergodicity parameter in the quiescent MSD and of the value for a. 
The agreement for long times is hence true without any fit parameter. As noted above, this 
agreement of the simulation data with Eq. (13-7bl) means that the long time value of the FDR 
takes the value of « \ and that the last term in Eq. (12-101 ) is indeed very small at long times. 
Future work has to show whether this is true in general. 
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§5. Summary 



The relation between mobility and diffusivity for shear melted glasses was presented. 
The derived relation followed without further approximation from previous results for cor- 
relation functions by taking the small wavevector limit. At short times, the Einstein relation 
holds. At long times, both diffusivity and mobility are finite and scale with shear rate, but 
the mobility is smaller than expected from the Einstein relation. In glasses, this scaling and 
the violation of the Einstein relation persist to arbitrary small shear rates, i.e., the limiting 
value of the FDR at 7 — > jumps at the glass transition from its nontrivial value X^ ~ ^ to 
its equilibrium value = 1. Within MCT-ITT, this jump has the same origin as the jump 
of the shear stress at the glass transition: The decay of the correlators on timescale I7I -1 in 
glasses. 

We find very good quantitative agreement to simulations. Future work has to address 
the only unknown term in the mobility, which is connected to the force-force correlation 
function. 
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Appendix A 

■ MCT-ITT equations for density correlators 



In this appendix, the MCT-ITT equations necessary to solve Eq. (13- lb are presented. In 
this paper, we will work on a schematic level, where the q-dependence in our equations for 
the correlators is dropped. The q-dependent treatment will be presented elsewhere. 8 ^ The 
schematic equation of motion for the normalized coherent transient correlator C coh (t, 0) = 
<P coh (t) reads 19 ' (this schematic model is called -model) 



= <P cotl (t) + r <J <P coll {t) + J dt'm cotl (j, t - t')$ con {t') 



m co/l (7,t) 



1 + {it/% 



cohy 



(vl + 2Ale)<P coh {t) + v c 2 {$ coh {t)) 



c /sF,coh 1 



(A- la) 
(A- lb) 



with initial decay rate r. We use the much studied values v% = 2, v\ = v^i^^/vr, — 1) 
and take m coh (0, t) in order to calculate quiescent (7 = 0) correlators. 34 ' In glassy states 
(e > 0), the long time decay of (p coh (t) from the plateau down to zero happens on timescale 
7 _1 , 19) i.e., without shear it stays on the plateau. For fluid states (e < 0), the correlator 
is analytic in shear rate and one observes a competition between structural relaxation on 
timescale r a and shear induced relaxation on timescale I7I 
strain jt, at which effects of shearing start to become important. We will use 7 



The parameter ^ c ° h sets the 

coh 

c 



1. 



Non-Equilibrium Einstein relation under shear 
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In MCT, the incoherent dynamics is coupled to the coherent one, in the sense that the 
coherent correlator enters the memory function of the equation for the incoherent correlator. 
The equation of motion for the incoherent density correlator under shear is known, 7 ) the 
schematic version for C(t, 0) = <P{t) reads 7 ' 

= (p(t) + r |*(t) + J dt'm(j, t - t')<P{t') | , (A-2a) 

771(7, t) contains the product of incoherent and coherent correlators, manifesting the cou- 
pling described above. This is why we introduce the coherent correlator here despite aiming 
only at tagged particle quantities. The effect of shearing is incorporated in a similar fash- 
ion as in (|A lb| ) and we again use m(0, t) in order to calculate quiescent quantities. 35) The 
additional parameter v s describes the coupling between the tagged particle and the bath par- 
ticles, i.e., different size ratios can in principle be mimicked. The standard MCT-analysis 
for Eqs. dA-2a| ) and dA-2b| ), see e.g. Ref. 30), yields the relation between the plateau values 
jcoh an( j j Q f Qcoh an( j Q reS p ec tively, 35) / = 1 tzsk- The incoherent dynamics is 

VsJ 

decoupled for values v s < v c s = l/f coh . For v\ and v% as chosen above, we have v c s m 3 
and use v s = 5 which is well above v c s . For v s > v c s , <P(t) has similai - properties as <& coh {t) 
for glassy and liquid states as described above. 7) Also 7 C = 1. 

The memory function 777,0(7, i n Eq. (13- lb finally is the q — > limit of the q-dependent 
analog of 777(7, t), which can be taken smoothly. We use 7770(7, t) = m {i->t) since we 
cannot perform the q — > limit in the schematic representation. 



Appendix B 

Normalized shear stress a 

(7) 

In the spirit of the F 12 -model, we approximate the s-dependent normalized shear mod- 
ulus by the transient correlator, 15 )' ) 

(ye^°a y ) ^ $co h{s) ^ (B-l) 

\O~xyO~xy) J 

where we account for the different plateau heights of the respective normalized functions by 
setting Goo/ f coh ~ f • cr follows then with Eq. (13-31 ) and t w — > 00. 
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